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Abstract 

We study heavy light mesons in a QCD inspired quark model with 
the Cornell potential— + br + c. Here we consider the linear term 
br as the parent and —^§f- + c i.e. the Coloumbic part as the pertur- 
bation. The linear parent leads to Airy function as the unperturbed 
wavefunction.We then use the Dalgarno method of perturbation the- 
ory to obtain the total wavefunction corrected upto first order with 
Coulombic piece as the perturbation. With these wavefunctions , we 
study the Isgur-Wise function and calculate its slope and curvature. 
Keywords: Dalgarno method,Isgur-Wise function, slope and curva- 
ture. 

PACS Nos. 12.39.-x ; 12.39. Jh ; 12.39.Pn 

1 Introduction 

Considerable efforts have been made in understanding the physics of hadrons 
containing at least one heavy quark since long [1-9]. It is well known that the 
heavy quark symmetry in the heavy quark limit leads to a single form factor 
called as Isgur-Wise (I- W) function which can describe the heavy quark bile- 
near current matrix elements of weak decay.The basic ingradient of the I-W 
function is the hadronic wave function, the determination of which becomes 
such a crucial factor. The potential models for this purpose is quite helpful 
as they contain more input parameters and hence has its firm basis. 
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Under such circumstances the I-W function has been investigated [3-9] 
with considerable success of vahd degrees in different models. In the poten- 
tial models 'Cornell potential' is found to be more useful than the others. It 
leaves two options of choosing the parent (1) the Coulombic part — and 
(2) the linear potential br. 

The slope and curvature of I-W function with the coulombic potential as 
the parent has already been reported for different heavy-light flavour mesons 
[10-15], which however had certain limitations. In coulombic potential as par- 
ent and linear as perturbation, the value of slope ( p^) and curvature (C) 
was found to be too large in MS -scheme. Imposing V-scheme [14-18], with 
larger as the values were found to be improved [13] but still larger than the 
expectations. As an alternate approach, in the present work we choose linear 
term '6r' as the parent and Coulombic piece as the perturbation in finding 
the wave function. 

As usual,two body Schrodinger equation is used with the ground state 
Fermi-Breit Hamiltonian in the absence of contact term and with Coulom- 
bic perturbation ,the wavefunction corrected upto first order is obtained by 
using the Dalgarno method [lOl [19]. The relavistic effect is incorporated 
by using standard Dirac modification [20l El] in a parameter free way. These 
wavefunctions are used in the calculation of slope and curvature of Isgur - 
Wise function. 

The rest of the paper is organised as follows : section 2 contains the for- 
malism , section 3 the result and section 4 the conclusion and discussion. 



2 Formalism 

2.1 The wavefunction 

We start with the ground state (/ = 0) spin independent Fermi-Breit Hamil- 
tonian without the contact term given by [10, ITT] : 



H 



2/i 



4a, 



3r 



+ 6r + c 



(1) 
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so that 

f^' = -^ + c (2) 
can be treated as perturbation to the unperturbed Hamiltonian : 



= + (3) 

In eq.(l), the strong couphng constant connected to the potential is a function 
of the momentum as 



a, 



where nf is the number of flavour. The constant 'c' at its natural scale 
is taken to be 1 GeV [13] The two body non relativistic Shroedinger wave 
equation can be recasted as 

H\ilj>={Ho + H')\ij>= E\ij> (5) 

The unperturbed wave function coresponding to Hq are the Airy functions 
which after normalization can be written as : 

i^^^Hr) = ^M{2f,by^ + pon) (6) 

where pon are the zeros of the Airy function (pon) = 0, n = 1,2,3.. 
represent the principal quantum no. (of course for the ground state n=l) and 
N is the normalisation constant. 

The pon are given as [221 ES] : 



.37r(4n-l),2 
Pon = -[^^ '-]^ (7) 

The first order corection to wave function -0^^-' and energy W^^-* are re- 
spectively given by 

Ho^l^^ + = W^^' + W^ij' (8) 
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where W° is the unperturbed energy given as 



and 



< = En 



Jo 



POn 



^(0) 



dr 



(9) 



(10) 



Since we consider the ground state (n = 1 ) , so we drop the 'n' from W^,Wji^\ 
ipl^'^ and -0^. The first order correction is : 



(r) 



4a^ Mo 



(11) 



As Airy function Ai(r) involve infinite series in r, so in calculating the coef- 
ficients ao,ai and a2 we have considered upto order and are given by : 



ao 



0.8808 (bfiY 



(32 



+ 



AW^ X 0.21005 



ai 



- c) ^ - c) 3a, {E - c) 
6ao 4 X X 0.8808 X (6//)^ 0.6535 x 



+ 



^2 



3a;s {E — c) 

_ AfxW^ X 0.1183 
3a. 



(E-c) 



(12) 

(13) 
(14) 



The total wave function corrected upto first order with normalisation is 

i^couiir) = ^/^(°) (0 + ^^'UO (15) 



2v^ 



y4i((2//fe)a + poi) /oo 



+ fli + 



(16) 



where A^i is the normalisation constant for the total wave function i/jcoui [r) 
with subscript 'coul' means coulombic potential as perturbation. 

The relativistic version of eq (17)is obtained by multiplying it with (^ ^J ^ ) 

o-Bohr depends on as : 



4 



a-Bohr — -, (17) 



and 



^ = l-\/l-(^) (18) 



Thus,relativistic wavefunction is : 

Aei (r) = ipcoui (r) (^—] (19) 

\CtBohr/ 

2.2 Isgur-Wise function 

The Isgur-Wise function is written as [H |2] : 

= l-p^y-l) + Ciy-lf + ... (20) 

where 

y = v,.v'^ (21) 

and and v'^ being the four velocity of the heavy meson before and after 
the decay. The quantity is the slope of I-W function aX y = 1 and known 
as charge radius : 

/ - 1 1.- (22) 

The second order derivative is the curvature of the I-W function known as 
convexity parameter : 

For the heavy-light flavor mesons the I-W function can also be written as 

mm ■■ 



47rr^ lip {r)f cos prdr 



(24) 



where 



= 2^' {y - 1) (25) 

the wavefunction eq.(19) with relativistic effect is used in the calculation 
of Hy) given by eq.(24). 



3 Calculation and Results 

We have calculated the values of charge radius and convexity parameter of 
the I-W function given by eq.(20) for two set of coupling constants both in 
WS and V -scheme [14-17]. 

Regarding the use of the above mentioned schemes [TT|, [T2] we note that 
with Uf = 4 and nj = 5 and fixing Aqcd = 0.2166*61^ [24]the coresponding 
value of cejjs at the scale of 1.5 GeV and 8 GeV are respectively 0.39 and 0.22 
[21] . The respective change of ajfg to (^) in the V scheme [H [151 M [H] 
for three different choices of scale JI are calculated [12] and shown in the ta- 
ble 1 here. Although there is no fundamental reason for the choice ,we have 
chosen the two renormahzation schemes {MS and l^-schemes) to fecilitate 
the comparision of our result with the previous work[12,13]. Also we use the 
same model parameter b = O.lSSGeV^ from charmonium spectroscopy[25,26] 

For these calculations, we have used the expressions for E, ,ao,ai given 
by equations (10), (12), (13), (14) respectively. These are shown in the table 2 
and table 3. The result of and c in the present work is shown in table 4. 
We also compare the present result with that of previous work with linear 
as the perturbation [13] in V^-scheme which was an improvement over MS 
-scheme and is shown in table 4. 

In table 4, we give a list of predictions of and C in different theoretical 
models. 

In evaluating the various integrations, we use numerical method of inte- 
gration in mathmatica software. 
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Table 1: The value of {^) for different choices of . 



choices 


7i=i 




"1 


^Msi'^b) = 0.22, n/ = 5 


0.259 


0.261 


0.258 


ajjs {rric) — 0.39, n/ = 4 


0.693 


0.651 


0.604 



Table 2: The values of and E in GeV 



Mesons 


E 




MS scheme 


V scheme 


D 


0.3898 


0.0467 


0.08314 


Ds 


0.4291 


0.5137 


0.0915 


B 


0.4072 


0.02742 


0.0327 


Bs 


0.4553 


0.0308 


0.0366 



Table 3: List of ao,ai and 02 



Mesons 


ao 


ai{GeV) 




V scheme 


MS scheme 


V scheme 


MS scheme 


V scheme 


MS scheme 


D 


0.2143 


0.1943 


-0.006138 


-0.007877 


0.00293 


0.002933 


Ds 


0.238 


0.21387 


-0.00916 


-0.01257 


0.0043 


0.0043036 


B 


0.2245 


0.2029 


-0.00749 


-0.0099 


0.00349 


0.00348 


B, 


0.251 


0.2269 


-0.0111 


-0.01601 


0.005116 


0.005 17 
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Table 4: Values of and C in our work and its comparision to other work 



Our work 


Scheme 


Mesons 




C 


MS'-scheme 


D 


0.7936 


0.0008 




Ds 


1.186 


0.002 




B 


0.89 


0.0004 




Bs 


1.41 


0.0012 


V- scheme 


D 


0.896 


0.00306 




Ds 


1.352 


0.0077 




B 


0.912 


0.0007 




Bs 


1.421 


0.00155 


Other work 


Previous work |1 2 1 113] 


D 


1.136 


5.377 




Ds 


1.083 


3.583 




B 


128.28 


5212 




Bs 


112.759 


4841 


Le Youanc et al [27j 




> 0.75 




Le Youanc et al [28] 




> 0.75 


> 0.47 


Rosner ^ 




1.66 


2.76 


Manuel [301 EI] 




0.98 


0.98 


Pole Ansatz [32] 




1.42 


2.71 


Ebert et al [36] 




1.04 


1.36 


Simple Quark Model [3] 




1 


1.11 


Skryme Model [35] 




1.3 


0.85 


QCD Sum Rule [H] 




0.65 


0.47 


Relativistic Three Quark Model [1] 




1.35 


1.75 


Neubert [33] 




0.82±0.09 
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4 Discussion and Conclusion 



Our calculated values of slope of I-W function in this work are found to be in 
good agreement with the other theoretical results (table 4). The lattice QCD 
evaluation of = 0.83llij;22 for B meson[37] and the experimental values of 
D meson pi = 0.76 ± 0.16 ± 0.08 [38] and pi = 0.69 ± 0.14 [39] are also in 
good agreement with our calculated results. However, the values of C for each 
meson are found to be smaller in comparison to other theoritical values. The 
reason may be presumably due to the cut off of the infinite series of Ai(z) 
upto O(r^) as noted earlier and still such small values can be considered as 
a success particularly for the B sector mesons as these values were very large 
in case of coulombic potential as parent [11- 14]. 

This study of the Isgur-Wise function with Coulombic part as perturba- 
tion shows a different picture as compared to the earlier work [111 [I2l [13] . 
With linear part as perturbation, the slope and curvature decrease with the 
increase of as] while in this work , we have observed a reverse effect. Fur- 
ther, this analysis shows a great reduction in the values of and C for all the 
mesons as compared to the previous work with linear part as perturbation. 
Let us conclude the section with a few comments. 

The strong coupling constant entering the coulombic potential is a function 
of the momentum in full QCD. But in potential model, it is nothing but a 
mere parameter. Here we have used the strong coupling constant in the MS 
and V^-scheme to fecilitate a proper comparision with the previous work with 
linear part as perturbation[12,13]. 

However, instead of using a particular renormalization scheme we could as 
well have considered the strong coupling constant merely as a free parameter 
in the potential model to be fitted from data. Such a possibility is currently 
under study. 
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